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Abstract 

The weak limit of the sequential empirical copula process is obtained and the 
asymptotic validity of a resampling scheme based on multipliers is established under 
the nonrestrictive smoothness conditions considered by Segers (2012). The empiri- 



cal process under consideration differs from the sequential process initially studied 



by Riischendorf (1976) which cannot be expressed in terms of the empirical copula. 



The obtained theoretical results are used to derive tests for detecting distributional 
changes in a sequence of independent continuous observations. The tests are based 
on the comparison of the empirical copula of the k first and the n — k last obser- 
vations for all possible values of k. The finite-sample performance of the resulting 
testing procedures is studied through large-scale Monte Carlo experiments. In the 
case of distributional changes due to a change in the copula only, the derived tests 
appear to be substantially more powerful than similar tests recently considered in 
the literature based on the sequential process initially studied by Riischendorf. 

Keywords: breakpoint detection; empirical copula; multiplier central limit theorem; 
multivariate independent observations; partial-sum process; ranks. 



1 Introduction 



Let X be a <i-dimensional random vector with continuous marginal cumulative distribu- 
tion functions (c.d.f.s) Fi,... , F d . It is then well-known from the work of Sklar (1959) 
that the c.d.f. F of X can be written in a unique way as 



F(x) = C{F 1 (x 1 ),...,F d (x d )}, xe 



1 



where the function C : [0, Vf — > [0, 1] is a copula and can be regarded as capturing the 
dependence between the components of X. 

Assume that C and F 1 , . . . , F d are unknown, and let Xi = (X a , . . . , X id ), i G {1, . . . , n}, 
be a random sample from F. For any j G {1, . . . , d}, denote by R%^ n the rank of X^ 
among Xij, . . . , X n j and let Uij, n = Rij,nl (n + 1). The random vectors 



Ui, n = (Ua jn , ... , U idi n), i G {1, . . . , n}, 
are frequently referred to as pseudo-observations from the copula C. 



(2) 



A natural nonparametric estimator of C is then the empirical copula of X\, . . . , X n 



( Riischendorf 1976 Deheuvels, 1979, 1981), which is frequently defined as the empirical 
c.d.f. computed from the pseudo-observations, i.e., 

1 n 

c r n («) = -2^i(t/i, n <«) J «G[o,if. 



i=l 



Note that the quantities Uij, n can equivalently be rewritten as Uij, n = nF n j(Xij)/(n + 1), 
where F n j is the empirical c.d.f. computed from X±j, . . . ,X n j, and where the scaling 
factor n/ (n + 1) is classically introduced to avoid problems at the boundary of [0, l] d . 

The object of interest of this paper is the sequential empirical copula process defined 

by 

[ns] 



CJs,u) 



V i=l 



< U) 



C(u) 



(s,u) G [0,1] 



d+l 



(3) 



where, for any y > 0, \_y\ is the greatest integer smaller or equal than y, and, for any 
s G [0, 1], the Ui^ ns \ are pseudo-observations computed analogously to (2) with n replaced 
by \ns\ . The latter process can be rewritten in terms of the empirical copula C[ ns j of 

Ax, . . • , Xyn^ as 



CJs, u) 



u 



C(u)}, 



(s,u) G [0,1] 



d+l 



where A n (s) = \ns\/n and with the convention that C (u) = for all u G [0, l] d . 
For s — 1, one recovers the standard empirical copula process which has been extensively 



g- 


Ruschendorf[ 


1976 


Ganssler and Stute 


1987 


Fermanian 


van der Vaart and Wellner 


, 2007; 


Sogers. 


20121 


Biicher and 



Volgushev, 2011). Note that C n as defined in ^ differs from the sequential process 



initially studied by | Ruschendorf| ( |1976 [) and defined by 

C(u)}, 



C*(a,«) 



LrtsJ 
i=l 



< U) 



(s,u) G [0,1] 



d+l 



(4) 



Unlike C n , cannot be rewritten in terms of the empirical copula unless s — 1. The 



latter 


was further studied by 


Biicher and Volgushev 


(2011) and was used in 


Remillard 


(2010 


), 


Biicher and Ruppert 


( 


2013 


), 


Wied et al. 


(2013) 


and 


van Kampen and Wied ( 


2012) 



to derive tests for change-point detection. 
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The first aim of this paper is to study the asymptotics of the sequential empirical 
copula process C n defined in ^ under the nonrestrictive smoothness conditions consid- 
ered in Segers (2012), and to show, in this more general framework, the validity of the 



multiplier bootstrap initially proposed by Remillard and Scaillet (2009). The second ob- 



jective of this work is to apply the obtained theoretical results to the derivation of tests 
for change-point detection for serially independent observations particularly sensitive to 
changes in the copula. 

Given a sequence Xi, . . . , X n of multivariate independent observations, change-point 
detection classically aims at testing 



Hq : 3 F such that X\, . . . , X n have c.d.f. F 



(5) 



against alternatives involving the nonconstancy of the c.d.f. Under the assumption that 
Xi, . . . , X n have continuous marginal c.d.f.s, from ([I]), we can write H as H 0>m n H c , 
where 



and 



ifo,m : 3 Fx, . . . , Fa such that Xi, . . . , X n have marginal c.d.f.s Fi, 



Hq c : 3 C such that Xi, . . . , X n have copula C. 



(6) 



(7) 



Following Csorgo and Horvath (1997, Section 2.6) (see also Bai, 1994), a broad class 



of nonparametric tests for testing H can be derived from the process 

B*(s,x) = V^K{s) {1 - X n (s)} {F [nsi (x) - FZ_ [nsi (x)} , (s,x) e [0, 1] x R d , 
where 



F,(x) 



1 k 

i=l 



< x) 



and 



i - 

zrr. E '« 



n — k 



< x) 



x e 



i=k+l 



are the empirical c.d.f.s computed from X\, . . . , X/, and . . . , X n , respectively, with 

the convention that Fq(x) = and Fq (x) = for all x e R d . A typical test statistic is 
then 

k(n-k) ^| _ _ ,™ # . 



T 



max 

Kfc<n-1 



n 



3/2 



sup \F k (x) - F, 



n—k 



X 



sup sup |^ „ 

se[o,i] £ceM d 



s,x) 



which was for instance used in Gombay and Horvath (1999). While D* appears as 



a natural process for testing H , typical tests derived from D* (such as those based 
on T n ) appear to have little power, in the case of moderate sample size, for detecting 
distributional changes due to a change in the copula only, i.e., when H ^ m D (->Hq iC ) holds. 



Empirical evidence of the latter fact can be found in Holmes et al. (2013, Section 4) 
and will be additionally given in Section |4j For that reason, nonparametric tests for 
change-point detection particularly sensitive to changes in the dependence structure are 
of strong practical interest. 

Several such tests were proposed in the literature. Tests based on Kendall's tau were 



investigated by Gombay and Horvath (1999), Gombay and Horvath (2002), Quessy et al. 
d2013| and pehling et al.| Q2013| ). Tests based on a variant of process (|8j) related to 



the process defined in Q were considered in Remillard (I2010I), Biicher and Ruppert 
(]2013[), |van Kampen and Wied| (|2012[) and|Wied et al.|(|2013). While some of these tests 



were designed to deal with serially dependent data, to the best of our understanding, 
none of these tests can be used to directly test H Q c defined in ^ unless H 0jTn defined 
in ^ is assumed to hold. In other words, when using these tests, rejection of H should 
not be interpreted as rejection of H Q c unless one is ready to assume that H rn holds. 
This point will be discussed further in Section |5j 

The second aim of this work is thus to apply the theoretical results obtained for 
the process C n defined in ^ to derive procedures for testing H based on yet another 
variant of process ([8]). As we shall see in Section 3.2, there are theoretical reasons to 
expect that, at least in some scenarios under H Qtm , the proposed tests are more powerful 



than the similar tests considered in Remillard (2010) and Biicher and Ruppert (2013). 



The numerical experiments of Section [4] empirically confirm these theoretical hints and 
reveal that, in finite-samples, the gain in power can be substantial. In a related manner, 
Corollary [4] in Section 3.2 suggests that, under H 0t7n , should the marginal c.d.f.s Fx, . . . , Fd 
be known, it is actually preferable to ignore that information, at least in some scenarios, 



thereby providing yet another illustration of a phenomenon put forward in Genest and 



Segers (2010) and, more recently, in Biicher (2012). 



The paper is organized as follows. In the second section, the weak limit of the sequen- 
tial empirical copula process C n is obtained and the validity of a multiplier bootstrap 
for C n is established under smoothness conditions a la Segers (2012). The third section 



is devoted to a detailed description of the tests for change-point detection based on the 
results of Section [2] as well as to a theoretical comparison with the similar tests considered 
Remillard ( 2010[ ) and Biicher and Ruppert (2013). The fourth section partially reports 



in 



the results of large-scale Monte Carlo experiments comparing the derived tests with those 



of Biicher and Ruppert (2013) and a test based on Kendall's tau considered in Quessy 



et al. (2013), among others. All the proofs are relegated to the appendices. 



In the rest of the paper, the arrow 



' denotes weak convergence in the sense of 
Definition 1.3.3 in van der Vaart and Wellner (2000), and £°°([0, l] d+1 ) represents the 
space of all bounded real- valued functions on [0, l] d+1 equipped with the uniform metric. 

Note finally that the code of all the tests studied in this work will be released as an 
R package whose tentative name is npcp. 



2 The sequential empirical copula process 

Let Xi, . . . , X n be a random sample from c.d.f. F and recall that Fx, . . . ,Fd are the 
continuous univariate margins of F. Then, let Ui = (Ua, . . . , Uid), i £ {1, . . . , n}, be 
the unobservable random sample obtained from Xi, . . . , X n by the probability integral 
transforms Z7 i3 - = Fj(Xij), i £ {1, . . . , n}, j £ {1, . . . , d}. It follows that Ux, . . . , U n is 
a marginally uniform random sample from the unknown c.d.f. C . The corresponding 
sequential empirical process is then defined as 

\ns\ 

Z n (s, u) = — $"{1(17; < u) - C(u)}, (s, u) £ [0, l] d+1 . (9) 



4 



According to Theorem 2.12.1 of |van der Vaart and Wellner| ( |2000| , Z ra ~* Z c in £°°([0, l] d+1 ) 
where Zc is a tight centered Gaussian process with covariance function 



cov{Z c (s, u),Z c (t, v)} — (s A t){C(u A v) - C(u)C(v)} 



(10) 



known C-Kiefer-Miiller process. 

To state the weak limit of the sequential empirical copula process C n defined in ([3]), 
we need the partial derivatives of the copula C. For any j G {1, . . . ,d}, let C»\u) be 
the partial derivative of C with respect to its jth argument at u, i.e., 



C [j] (u) 



lim 

h->0 
u,-+hG[0,l] 



C(ui, . . . ,ttj_i,Mj + h,u j+1 , . . . ,u d ) - C(u) 
h 



u e [o, l]' 



If C^(u) exists and is continuous on [0, l] d for all j G {1, . . . , off, we know from the work 



of Ganssler and Stute (1987), Fermanian et al. (2004), Tsukahara (2005) or van der Vaart 



and Wellner (2007) that the empirical copula process y/n(C n — C) = C n (l, •) converges 
weakly in £°°([0, l] d ) to the tight centered Gaussian process 



C c (l, u) = Z c (l, u)-J2 C^(u)Z c (l, ««), u e [0, 1]' 

3=1 



'IF 



where, for any j G {1, . . . ,d} and any u G [0, l] d , is the vector of [0, l] d defined by 



u 



U) 



Uj if i = j and 1 otherwise, and Zc(l, ■) is a C-Brownian bridge. 



For many copula families however, the partial derivatives j G {l,...,d}, fail 
to be continuous on the whole of [0, l] d . To deal with such situations, Segers (2012) 
considered the following condition: 

(CI) for any j G {1, . . . , d}, exists and is continuous on the set Vj = {u G [0, l] d : 
< Uj < 1}. 



Under Condition (CI), for any j G { 1 , . . . , d], Segers (2012) extended the domain of 
to the whole of [0, l] d by setting 



C [j] ( 



C(ui, . . . , uj-i, h, u j+1} ...,u d ) 
h 



Hue [0, l] d , Uj = 0, 



u 



lim sup 

hmsup C{u)-C{u^...^l-Ku 3+U ..^u d ) ^ . f u e [Q) = ^ 

Ho " 



which ensures that the process Cc(l, •) defined in (11) is well-defined on the whole of 



[0, l] d , and showed the weak convergence of the empirical copula process \fn{C n — C) = 
C„(l, ■) to Cc(l, ■) in £°°([0, l] d ). Condition (CI) is nonrestrictive in the sense that it is 



necessary for Cc(l, ■) to have continuous sample paths and was verified in Segers (2012) 
for many popular copula families. 

A stronger result can be obtained if the following additional smoothness condition is 
assumed: 
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(C2) for any i,j G {1, . . . ,d}, the second-order partial derivative C' 1 '^ = d 2 C /(duiduj) 
exists, is continuous on the set V{ fl Vj — {u G [0, l] d : Ui, Uj G (0, 1)}, and 

sup max{iij(l — Ui),Uj(l — Uj)}\C^^\u)\ < oo. 



Under Conditions (CI) and (C2), Segers (2012, Proposition 4.2) showed that, with prob- 
ability one, 



sup 

■ue[o,i] a 



C„(l,u) - \ Z n (l,u) -J2C [j] W n (l,u^) 



0(n- 1/4 (logn) 1/2 (loglogn) 1/4 ), 

(12) 



thereby recovering the result initially stated by Stute (1984) under very restrictive con- 
ditions on C . 



The result stated in (12) is at the root of the weak convergence of the sequential 



empirical copula process C n defined in fl3J). The proof of the following proposition is 
given in Appendix |A| 

Proposition 1. Under Conditions (CI) and (C2), C n -w Cc in £°°([0, l] d+1 ), where 



C c (s, u) = Z c (s, u)-J2 C bl Wc(s, u 



(s,u) e [0,1] 



d+l 



(13) 



The computation of approximate p-values for tests involving the sequential empirical 
copula process will typically be based on some resampling scheme. In the nonsequen- 



tial setting, Biicher and Dette (2010) compared the finite-sample behavior of the various 



resampling techniques proposed in the literature and concluded that the multiplier pro- 



cedure of Remillard and Scaillet (2009) has, overall, the best finite-sample properties. 



This technique was revisited theoretically by Segers (2012) who showed its asymptotic 



validity under Condition (CI). In the remainder of this section, we show the asymptotic 
validity of the multiplier approach in the sequential setting under consideration. 

Let M be a large integer and let £ f , i G {1, . . . , n}, m G {1, . . . , M}, be i.i.d. random 
variables independent of the data Xi, . . . , X n . We consider two possibilities for the £f m ^: 



(All) Q has mean 0, variance 1 and satisfies / ° O {P(|6j m ''| > x^^dx < oo; 



Mi 



(A42) £^ is a Rademacher random variable, i.e., P(£ 



(m) 



(m) 



■1) = 1/2. 



Clearly, a result that holds under will also hold under (A42). 

Next, for any m G {1, . . . , M}, let 



s,u 



|nsj 

z= ct ] im < «) - c(u)}, ( s , u) g [o, i] 



d+l 



(14) 
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A consequence of the sequential extension of the multiplier central limit theorem proved 
in 



Holmes et al. (2013, Theorem 1) is that 



(M) 
C 



in {£°°([0, 1] )} M+1 , where Zc, a C-Kiefer-Miiller process, is the weak limit of the 
sequential empirical process Z n defined in (|9), and Z^, . . . , Zj^ are independent copies 
of Zc- For large n, Z„ , . . . , Zn could therefore be regarded as "almost" independent 
copies of Z„. Unfortunately, the zi" 1 '' cannot be computed because the random sample 
Ui, . . . , U n is unobservable and C is unknown. 

We consider two versions of the Zl m) that can be computed depending on how the Ui 
and C are estimated. For any s G [0, 1] and u G [0, l] d , let 



Zi m) ( 



s, u) 



[ns\ [ns\ 

-^£fc (m) {l(0U».j <^ 
* i=i * i=i 

where the E^Lnsj are defined as in (2) with n replaced by \ns\, = [wsJ -1 Yli=i 
and £q = by convention, and let 



|nsj 

V i=l 



(16) 



To define "almost" independent copies of C n for large n in the spirit of Remillard and 



Scaillet (2009), we also need to estimate the partial derivatives appearing in (13). To 



do so, we consider the estimator C$ defined in 
1 



Kojadinovic et al. (2011) by 



u 



— {C n (ui, . . . , Uj-i, uj n , Uj+x, ...,u d ) 

-C n {u l ,... 1 u j ^u~ n1 u j+l} ...,u d )} , ue[0,l} d , (17) 



where u^ n = (uj + n~ 1 / 2 ) A 1, and u7 n = (uj — n" 1 / 2 ) V . This estimator differs slightly 
from the one initially proposed in Remillard and Scaillet (2009). It has the advantage of 
converging in probability to uniformly over [0, l] d if happens to be continuous 
on [0, l] d instead of only satisfying Condition (CI). 

The following result, proved in Appendix [A] complements Lemma 2 of Kojadinovic| 



et al. (2011). 



Proposition 2. Let < a < b < 1 and assume that Conditions (CI) and (C2) hold. 
Then, for any j G {1, . . . , d}, 



sup \CH ] (u) - C®(u)\ = O (n- 1/4 (logn) 1/2 (loglogn) 1/4 ) 

u.£[0,l] d 



almost surely. 
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We can now define empirical processes that can be fully computed and that, under 
appropriate conditions, can be regarded as "almost" independent copies of C n for large n. 
For any m G {1, . . . , M}, s G [0, 1] and u G [0, let 



Ci m) ( 



Zl m) ( 



and 



(19) 



The following result, proved in Appendix [Bl can be seen as an extension of Proposi- 



tion 3.2 of Segers (2012) to the sequential setting. 



are 



Theorem 1. Assume that Conditions (CI) and (C2) hold, 
(i) Under (Ml), 

(C n ,C^\...^)^(c c X^\---X { c M) ,, 

in {£°°([0,l]( d+1 ))} M+1 ; where C c is the weak limit of C n and Cff, . . . , Cjf } 
independent copies Cc- 

(ii) If (M.2) is assumed instead of (Ml), 

(c n , c«, . . . , cf, c«, . . . , c[ M >) - (c c , cg } , . . . , cjf \ cg\ ... , cg M) ) , 

m {£°°([0,l]( d+1 ))} 2M+1 . 



3 Tests for a change in the dependence structure 

The setting is as follows. Given a sequence of independent continuous ci-dimensional 
random vectors Xj, . . . , X n , the objective is to test H defined in ^ against alternatives 
involving the nonconstancy of the c.d.f. The aim being to derive tests particularly sen- 
sitive to changes in the copula of the sequence, the behavior of the proposed tests will 
be theoretically investigated under the alternative hypothesis of constant marginal c.d.f.s 
and a single change-point in the copula, i.e., under H 0>m fl Hi tC , where H 0tm is defined 
in ^ and Hi c is given by: 

Hi c : 3 distinct C\ and C*2, and k* G {1, . . . , n — 1} such that 

X±, . . . , Xk* have copula C\ and Xk*+i, • • • , X n have copula Ci- (20) 



After describing the proposed tests, we will compare them with the similar tests 



considered in Remillard (2010, Section 5.2) and Bucher and Ruppert (2013 Section 3.2) 



and provide theoretical hints suggesting that, under Ho iTn , the former should be more 
powerful than the latter at least in some scenarios. 



S 



3.1 Tests based on the empirical copula 

As explained in the introduction, the aim is to adapt to the current copula context the 



approach studied in detail in Csorgo and Horvath (1997, Section 2.6) based on the process 
defined in ([8]). The idea consists of comparing, for all k G {1, ... ,n — 1}, the empirical 
copula of Xi, . . . , Xk with the empirical copula of Xk+i, • • • , X n , respectively defined by 



CM 



1 k 



< u) 



i=i 



1 

and C* n _ k {u) = — ^tn-u < u), u G [0, 1]' 



i=k+l 



X r , 



where U k , ln _ k ,...,U* n _ k are the pseudo-observations computed from Xk+i, 
similarly to ([2]). Analogous to ([8]), we define the process 

B n (s,u) = V^K(s) {1 - X n (s)} {C [nsj (u) - C*_ lns] (u)}, (s,u) G [0, l] d+ \ (21) 

where A n (s) = \ns\/n and with the convention that Co (it) = and Cq(u) = for all 
it G [0, l] d . 

The weak limit of D n is established in the following result proved in Appendix |Cj 

Proposition 3. Under Conditions (CI) and (CI), and H , D n ~* D c in £°°([0, l] d+1 ), 
where, for any (s,u) G [0, l] d+1 , H>c(s,u) = Cc(s,u) — sCc(l,it), with Cc defined 



in (13), and where the covariance function of "Be can be expressed in terms of the covari- 
ance function of Cc as 



cov{D c (s, u),B c (t, v)} = (s A t - st) cov{C c (l, u), C c (l, v)}. 



(22) 



To compute approximate p- values for tests based on D„, we adapt the multiplier 
approach described in the previous section to the current change-point setting. Analogous 
to Zn and defined in (15) and (16), respectively, for any m G {1, . . . , M}, s G [0, 1] 
and u G [0, l] d , let 



V„_[n«J — U ) ~ ^n-\ns\ 



* i=[ns\+l 
1 n 

* i=[ns\+l 

where £%-[ns} = ( n ~ [ns])' 1 YH=[ns\+i ^ and £o' (m) = b y convention, and let 



1 n 

K^\s,u) = -= ^ ] {m,n<u)-C n (u)}=Z^\l,u)-Z^(s,u) 



i=[nsj+l 

Then, for any s G [0, 1] and it G [0, l] d , let 



s, it 



£*,(m) ( 



3=1 



where C^t, is the version of the estimator of defined in (17) computed from 
X 

|nsj+i> • • • ) 3C n , and let 

d 

C*/ m \s, u) = Z*/ m \s, u)-J2 C)i\u)frf™\s, 

i=i 

Finally, for any m G {1, . . . , M}, s G [0, 1] and u G [0, l] d , we form 



„ { s,u) = {l-\ n (s)}C^\s,u)-\ n (s)C*/ m \s,u), 



and 



where and Cn' 1 ' are defined in (jlSh and (|19|), respectively. 



,(m) 



The following result, proved in Appendix [Cj is a consequence of Theorem [TJ 
Proposition 4. Assume that Conditions (CI) and (C2) hold. 

(i) Under (Ml) and H , 



in {£°°([0, l]( d+1 ))} M+1 ; where D c is the weak limit of B„ and o£ } , . . . , B[f ) 
independent copies B>c- 

(ii) If (M2) is assumed instead of (Ail), then, under Hq, 



are 



D D (m) D (1) 



The next result, also proved in Appendix [Cj will be useful to derive the asymptotic 
behavior of tests based on D n under H Q m PI Hi c , where H m and H\ )C are defined in (|6j 
and (20), respectively. 



Proposition 5. Assume that H 0rn r\Hi c holds with k* = \nt\ for some t G (0, 1). Then, 



(i) n-^ 2 B n A K t m£°°([0,l} d+1 ), where K t (s,u) = (s At)(l - s Vt){d(u) - C 2 {u)}, 

sup se[01] sup ue[0jl]d |M m) 
for all m G {1, . . . , M}. 



(ii) sup se[01] sup u6[01]d \B ( n n \s,u)\ = P (1) and sup s6[01] sup ue[0)1]d \B ( ™\s,u)\ = P (1) 



In the framework under consideration, a change in the c.d.f. of the sequence X\, . . . , X n 
can occur at any point k G {1, . . . , n — 1}. As classically done, a test for change-point 
detection can therefore be obtained by first defining a test statistic for any possible 
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change-point k G {1, . . . , n — 1}, and then by combining the resulting n — 1 statistics into 
a global statistic. 

Natural choices for the n— 1 change-point statistics are Cramer-von Mises, Kolmogorov- 



Smirnov and Kuiper statistics (see e.g. Horvath and Shao, 2007 Biicher and Ruppert 



2013). As the latter two statistics led to consistently less powerful tests in our Monte 
Carlo experiments, for the sake of brevity, we only describe the tests based on the n — 1 
Cramer-von Mises change-point statistics defined by 



S„ 



n,k 



[0,1]« 



k 



n 



u) \ dC n {u), k G {l,...,n- 1}. 



(23) 



A natural way of combining S^i, . . . , S n>n ^x consists of taking their maximum. This leads 
to the maximally selected Cramer-von Mises statistic 



S n 



max 

Kk<n- 



S n , k = sup / {D n (s,u)} 2 dC n (u). 
se[o,i] J[o,i] d 



(24) 



Note that we did also consider the arithmetic mean as combining function as in Holmes 



et al. (2013) but do not insist on this point further as the resulting tests did not appear 



more powerful than those based on S n in our Monte Carlo experiments. 

We studied two versions of the test based on S n depending on whether the approximate 
p- value is computed using the multiplier processes 6l" 1 ' ) or the processes lD>i m \ For any 
m G {1,...,M}, let 

£i m) =sup / (s,u)Y dC n (u) and = sup / 0, u)f dC n (u). 

se[o,i] J[o,i] d L J *e[o,i] ^[o,i] d 

The following result is then a corollary of Proposition H] and can be proved along the 



lines of Proposition 1 in Holmes et al. (2013). 



Corollary 1. Assume that Conditions (CI) and (CI) hold, 
(i) Under (Ml) and H , 

in [0, oo)( M+1 \ where 

5= sup / {D c (s,u)} 2 dC(u) 

ae[0,l] J[0,l} d 

is the weak limit of S n , and . . . , are independent copies of S. 
(ii) If (AA2) is assumed instead of (Ail), then, under H , 

(q A(l) A(M) r-(l) q(M)\ (q q(l) q(M) <j(l) q(M)\ 

m [0,oo)( 2M+1 ). 
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The previous corollary suggests interpreting the (resp. the Sn) under the null 
hypothesis as M "almost" independent copies of S n and thus computing an approximate 
p- value for S n as 

M M 



or as 



m=l 



(25) 



m=l 



The following result is a corollary of Proposition [5] and can be proved along the lines 
of Proposition 2 in Holmes et al. (2013). 

Corollary 2. Assume that H rn D ifi jC /ioWs wift k* = \nt\ for some t e (0, 1). JTien, 
{B n (s,u)} 2 



(i) s i — v 



-dCJu) 



[0,1]' 



n 



[ {^(^D^w + (l -t)c 2 (u)} 

J[0,l] d 



in £°°([0, 1}), which implies that S n — > +oo, 

(ii) for any me {1,..., M}, = P (1) and = O p (1). 

A consequence of the previous corollary is that, under i?o,m H i?i )C , the approximate 
p- values for 5* n defined in ( 25 ) will tend to zero in probability 

Finally, should H be rejected, and assuming that H 0t7n D H\ fi holds with k* = \nt\ 
for some t e (0, 1), it is frequently of interest in applications to estimate t. In the setting 
under consideration, an estimator of t can be defined as t n = n -1 argmax 1<A;<n _ 1 S nj k, 
where the S n ^ are defined in (23). 

Corollary 3. Assume that H 0tm PI Hi c holds with k* = \nt\ for some t e (0, 1). Then, 

J- P , 4. 

t n -»> *• 

The above result follows from Corollary [2] (i), the fact that, as a function on [0, 1], 
s i — y Jj Q ^ d {K t (s, u)} 2 d{tCi(u) + (1 — t)C 2 (u)} attains its unique maximum at s — t, and 



the argmax theorem (van der Vaart and Wellner, 2000, Corollary 3.2.3); see also Bucher 



and Ruppert (2013 Proof of Proposition 3) and Dehling et al. (2013 Proof of Theorem 



3.2 Comparison with tests related to the sequential empirical 
process of Riischendorf 

Similar tests to those derived in the previous subsection were proposed by Remillard] 
(2010, Section 5.2) and Bucher and Ruppert (2013, Section 3.2) using the empirical 
process 



s,u) = yfn~\ n (s){l-\ n (s)}{C lnsi (u)-C*_ [nsj (u)}, (s,u) e [0,1] 



d+l 



(26) 



where, for any k £ {1, . . . , n — 1}, and C*_ k denote respectively the empirical c.d.f.s 
of the k first and n — k last pseudo-observations Ui tn , . . . , U n ,n defined in (0), i.e., 

k 



C,(u) 



i=l 



< u) 



and C*_ k {u) = V l{U iin <u), uE [0, l] c 

n k 



=k+l 
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and with the convention that Cq(u) = and Cq(u) = 0. Starting from O^, these authors 
naturally suggested to consider maximally selected Kolmogorov-Smirnov, Kuiper and 
Cramer-von Mises test statistics. The latter is given by 



S 



R 



sup 

se[o,i] J [0,1 



/ <' 

J\o,i] d 



s, u 



)} 2 dC n 



u 



(27) 



and is the analogue of the statistic S n defined in (24) from the process D n defined in (21 ). 



Actually, all the test statistics considered in Remillard (2010, Section 5.2) and Bucher 



and Ruppert (2013, Section 3.2) are rank-based, and, consequently, can be rewritten in 



terms of the process D* defined in pi). For instance, in the case of one obtains 



S,, 



sup 

se[o,i] 



p / {B*(s,x)} 2 dF n 

,11 JR d 



X 



(2f 



where F n is the empirical c.d.f. computed from X\, . . . ,X n . It follows that the corre- 
sponding tests should not be regarded as tests for change-point detection focusing in 
particular on changes in the copula, but rather as general-purpose tests for detecting 
changes in the multivariate c.d.f. 

The finite-sample performance of the test based on S„ (seen as a functional of B#) 



was actually also investigated in Holmes et al. (2013, Section 4). The test was found 
to have good power for detecting changes in the marginal c.d.f.s. but did not appear to 
detect changes in the copula well. Similar empirical findings will be presented in the next 
section. 

Let us now assume that H m defined in ^ holds. In the remainder of this section, 
we provide theoretical hints suggesting that, at least in some situations, tests based on 
D n should be more powerful than analogous tests based on D^. The main ingredient of 



our argumentation is a covariance inequality obtained by Genest and Segers (2010). 



m 



From Theorem 5 of Bucher and Ruppert (2013), we have that, under Hq 
i°°([0, l] d+1 ), where B£(s, u) = Z c {s, u) - sZ c (l, u) for all (s, u) e [0, l] d+ \ and Z c is 



a C-Kiefer-Miiller process whose covariance function is given in (10). Simple calculations 
show that the covariance function of can be expressed in terms of the covariance 
function of the C-Brownian bridge Zc(l, •) as 



cov{I 



(t, «)} = (s A t - st) cov{Z c (l, u), Z c (l, «)}. 



Combining the above expression, (22) and Proposition 1 of Genest and Segers (2010), we 
obtain the following result. 

Corollary 4. Assume that Condition (CI ) holds and that C is a bivariate copula that is 
left-tail decreasing in both arguments. Then, under Hq, for any (s,u), (t,v) 6 [0, 1] 3 7 



cov{Bc(s,u),B c (t,v)} < cov{B^(s,u),B^t,v)}. 



(29) 



From Nelsen (2006, Section 5.2.2 and Exercise 5.35), a bivariate copula C is left-tail 



decreasing (LTD) in both arguments if, for any < u < vl < 1 and < v < v' < 1, 



C(u,v)/(uv) > C(u',v')/(u'v'). As shown by Garralda-Guillem (2000), extreme- value 
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copulas are LTD in both arguments but so are the most popular bivariate copulas with 
positive dependence such as the Clayton, Frank, normal, t and Plackett. 



To illustrate a possible consequence of (29) and provide some intuition, assume that 
d = 2, that Condition (CI) holds and that C is LTD in both arguments. Furthermore, 
suppose that test statistics are defined using the map T from £°°([0, l] 3 ) to R given by 
T(D) = L 1]3 D(s, u)dsdu, D G £°°([0, l] 3 ), and that large values of |TQD>£)| and |T(D„)| 
lead to the rejection of H . Then, from the continuous mapping theorem and (29), we 
obtain that, under H , T(B*) ~* T(Dg) and T(B n ) ~* T(D C ), where T(Bg) and T(D C ) 
are centered normal random variables with var{T(Dc)} < var{T(D^)}. It follows that 
there exists xq > such that, for any x > xq, 



P{|T(D c )|>x}<P{|T(Dg)|>*}. 



(30) 



Now, assume that H 0in fl H l c holds with H l c defined in (20), k* 
t G (0, 1) and that the copulas C\ and Ci are such that J, Q 1 , 2 {Ci(u) — C2(u)}du ^ 0. 



\nt\ for some 



Then, from the proof of Proposition 3 of Biicher and Ruppert (2013) and from Proposi- 
tion [5] (i) in this work, both n~ 1//2 B^ and n _1 / 2 O n converge in probability in £°°([0, l] 3 ) 
to the constant function K t defined in Proposition [5] (i), and, by the continuous mapping 
theorem, both ra-^TQD*) and n" 1/2 T(D n ) converge in probability to T(K t ) 7^ 0. From 



the previous fact and (30), we can therefore expect that, as n increases, the asymptotic 
p-value for the test based on T(D n ) will tend to be smaller than the asymptotic p-value 
for the test based on T(D^). 



Unfortunately, at this point (see the discussion in Genest and Segers, 2010, Section 4), 
the previous argumentation cannot be extended to larger dimension or dependence con- 
ditions other than left-tail-decreasingness. Even for bivariate copulas that are LTD in 



both arguments, the case of the statistic S n defined in (24) and its analogue defined 



in (27) cannot be treated similarly, nor does it seem easy to deal with p- values computed 
approximately from multiplier replicates as in (25). Still, in our numerical experiments 



to be partially reported in the next section, the two versions of the test based on S n 



considered in Section 3.1 were found to be consistently more powerful than the analogue 
multiplier test based on S^. 

Let us end this subsection with a remark. Assume that H rn defined in ^ holds, 
and, for a moment, that we are in the unrealistic situation in which the continuous 
marginal c.d.f.s F\, . . . , Fa of the Xi are known. In this case, the sample Ui, . . . , U n 
with Ui = (Fi(Xn), . . . ,Fd(X i( i)) can be computed and it might seem intuitively better 
to use, instead of the process D n defined in (21), the process D* defined as the process 
but with the X t replaced by the E7j. Using rather classical arguments, it can 
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be verified that the difference between D* and defined in (26) converge uniformly in 



probability to zero, which implies that, under Hq, the two processes converge weakly to 
the same limit defined previously in this subsection. Corollary [4] then suggests that, 
at least in some scenarios, it is preferable to ignore the fact that the true marginal c.d.f.s 
are known and still use the process D n for testing. Besides Proposition 1 of |Genest and 



Segers (2010), yet another illustration of this phenomenon is given in Biicher (2012). 
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4 Monte Carlo experiments 



4.1 The setting 



Large-scale Monte Carlo experiments were carried out in order to study the finite-sample 
performance of various tests for change-point detection based on the empirical process D n 
defined in (21 . In the rest of the paper, the two versions of the test based on the statistic 



S n defined in (24) will be referred to as the test based on S n (resp. S n ) when approximate 
p- values are computed using the multiplier processes 3^ (resp. tl)n )■ As explained in 
the previous section, in addition to the tests based on the maximally selected Cramer-von 
Mises statistics S n and S n , we also considered tests based on Kolmogorov-Smirnov and 
Kuiper change-point statistics. Empirical results for these tests are not reported as they 
appear to be consistently less powerful than those based on S n and S n . 

The two aforementioned tests were also compared with the test based on the statistic 



defined in (27) derived from the process defined in (26). An approximate p- value 
for the latter was computed using the multiplier processes 



ivR,(m), 



s,u)E[0,l] d+1 ,me{l,...,M}, 



where Z 



(m) 



(16 


); see also 


Remillard 


(2010 



(2010, Section 5.2) and 



Biicher and 



Ruppert| p013l Section 3.2). 

As a second competitor to the tests based on S n and S n , we included a test based 



on differences of sample Kendall's taus. Such a test was initially studied in Gombay and 



Horvath 



Horvath 



(1999 



Section 3) in dimension two. Starting from the work of Gombay and 



(2002), multivariate extensions of this approach were recently investigated in 



Quessy et al. (2013). Given a ci-dimensional random vector X with continuous marginal 
c.d.f.s Fi, . . . ,F d and copula C, a natural multivariate extension of Kendall's tau is given 
by r = (2 d E{C(U)} 



l)/(2 



d-l 



1), where U = (F 1 (X 1 ), . . . , F d (X d )) (see e.g. Joe 



1997). One of the most powerful tests investigated in |Quessy et al. ( 2013[ ) is based on 
the statistic 

*(*-*)■- _* , (31) 



max 

2<k<n-2 



n 



3/2 



'n-kh 



where and r*_ fc are multivariate sample Kendall's taus computed from Xi, . . . , 
and Xk+i, • • • , X n , respectively. An approximate p- value for the latter test was obtained 



from multiplier replicates computed as explained in Gombay and Horvath (2002 page 
251) with kernel h defined as h(Xi,Xj) = < Xj) + > Xj) and up to the 

normalizing constant 2 d /(2 d ~ 1 — 1). 

To study the finite-sample performance of the tests based on S n , S n , and K n , data 
were first generated under H defined in ([5]), then under H , m nHi yCl where Ho jm is defined 
in (|6j and H\ jC is defined in (20), and finally under i?i >m fl Ho )C , where H c is defined 
in (fffi and 



H l rn : 3 distinct F^i, F 12 as well as F 2 , . . . , F d and k* £ {1, . . . , n — 1} 
such that Xi, . . . , X^* have marginal c.d.f.s F^i, F 2 , . . . , F d 
and X k * + i, . . . , X n have marginal c.d.f.s F 1>2 , F 2 , . . . , F d . 



(32) 
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Figure 1: The left (resp. right) plot shows 10,000 realizations of aC a ^ defined in (33) 
when C is the bivariate normal copula with parameter 0.97 and (a, (5) = (0.25, 1) (resp. 
(a,/3) = (0.5,1)). 



In terms of copula families, we considered the Clayton, Gumbel-Hougaard, Frank and 
exchangeable normal families. A nice overview of their properties is given for instance in 



McNeil et al. (2005, Chapter 5). To fix parameter values for data generation, we exploited 
the well-known one-to-one relationship between Kendall's tau and the parameter value 
in the bivariate versions of these models. To increase the diversity of shapes in dimension 
two, asymmetric versions of the previous copulas were defined using the simplest form 
of Khoudraji's device (Khoudraji, 1995 Genest et al. 1998 Liebscher 2008). Given an 



exchangeable bivariate copula C, Khoudraji's device can be used to define an asymmetric 
version of it as 



aC, 



a,/3{U,V) 



u 



U , V 



u,ve [0,1], 



(33) 



for arbitrary choices of a ^ (3, a, (3 6 (0, 1]. As an illustration, the left (resp. right) plot 
of Figure [T] shows 10,000 realizations of aC a ^ when C is the bivariate normal copula with 
parameter 0.97 and (a, (3) = (0.25,1) (resp. (a,/3) = (0.5,1)). 

Two choices were considered for the multipliers: Rademacher random variables and 
standard normal random variables. The two settings led to very similar results suggesting 
that Theorem [I] (ii) might actually hold under Condition (.Ml). The results reported in 
the forthcoming tables are those obtained with standard normal multipliers. 

For each data generating scenario, the rejection rates of the tests were estimated from 
1000 samples of size n G {50, 100, 200} and dimension d = 2 (and d = 3 when relevant) 



generated using the copula package ( |Hofert et al. , 2012) for the R statistical system ([Rj 
Development Core Team, 2012 ). All approximate p- values were computed from M = 1000 



multiplier realizations and the tests were carried out at the 5% level of significance. 
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Table 1: Percentage of rejection of i^o defined in ^ computed from 1000 random samples 

of size n e {50, 100, 200} generated under Hq, where C is either the d-dimensional Clayton 

(CI), the Gumbel-Hougaard (GH) or the normal (N) copula whose bivariate margins have 

a K endall's tau of r. 

CI GH N 



d n t S n S n K n S n S n Sff n n S n S n S^" 

1 50 0.00 6.2 4.0 4.6 5.8 5.4 2.9 4.5 3.6 7.3 3.4 4.8 5.1 

0.25 6.7 6.2 5.6 5.6 5.5 3.3 5.4 4.3 4.4 3.0 6.3 3.4 

0.50 5.6 7.9 6.0 3.1 4.4 3.3 4.6 3.6 4.4 5.3 4.9 2.9 

0.75 6.0 16.6 5.5 3.2 3.2 6.7 4.3 2.6 3.6 9.1 4.9 2.9 

100 0.00 4.9 3.5 5.3 5.2 5.2 4.1 5.5 5.8 4.3 2.8 5.5 3.8 

0.25 6.1 6.6 5.0 5.0 5.0 3.3 6.2 5.6 5.3 4.0 5.5 5.3 

0.50 4.4 9.3 5.9 4.1 3.7 2.8 5.7 3.6 3.1 3.4 5.3 2.8 

0.75 2.7 10.0 4.6 2.1 2.5 4.7 4.4 1.8 2.1 6.0 5.6 1.5 

200 0.00 4.0 3.5 5.2 4.4 4.3 4.0 5.2 4.4 5.4 4.9 4.3 5.3 

0.25 4.7 5.2 6.3 4.9 3.3 3.0 3.8 4.6 4.0 3.9 5.2 5.2 

0.50 5.1 8.5 4.9 5.0 3.2 2.3 4.5 4.4 4.0 4.7 4.8 4.2 

0.75 2.6 9.3 5.9 3.2 1.5 3.1 5.2 3.1 1.9 4.8 5.7 1.5 

3 50 0.00 4.3 1.5 3.0 2.6 4.2 2.1 3.6 3.4 5.5 2.8 3.4 4.4 

0.25 6.3 5.0 5.1 5.6 5.5 1.0 5.1 4.9 5.3 3.0 4.3 5.2 

0.50 7.2 8.1 5.9 6.2 2.7 0.9 5.7 4.1 3.0 2.2 4.6 3.6 

0.75 2.0 2.9 6.9 3.1 0.5 0.4 6.3 1.9 1.1 1.3 4.1 1.9 

100 0.00 4.5 3.4 4.5 4.3 4.5 2.8 4.6 3.7 4.5 2.7 3.9 5.0 

0.25 5.0 5.1 5.4 5.6 4.2 2.6 4.4 5.8 5.4 3.5 4.5 5.1 

0.50 5.7 7.6 6.3 5.0 3.3 1.3 5.0 4.6 3.2 3.1 3.9 3.6 

0.75 2.5 4.9 5.0 2.3 1.0 1.0 5.2 2.6 0.8 1.6 5.5 2.0 

200 0.00 3.3 2.5 4.3 4.1 3.5 3.2 4.3 4.1 4.8 4.0 4.7 5.2 

0.25 6.6 7.1 5.5 6.1 5.0 3.3 4.5 5.0 4.8 4.1 5.0 5.2 

0.50 6.0 9.2 4.5 5.6 3.0 2.4 5.9 5.3 4.8 4.3 4.8 5.5 

0.75 2.9 6.4 6.4 4.1 0.7 0.9 3.8 3.6 1.3 2.2 4.9 3.9 



4.2 Rejection percentages under Hq 

Table [I] reports the rejection percentages of Hq defined in ^ computed from 1000 
random samples generated under H , where the copula of F is either the Clayton, 
the Gumbel-Hougaard or the normal whose bivariate margins have a Kendall's tau 
of t G {0,0.25,0.5,0.75}. Since all the tests under consideration are rank-based, the 
marginal c.d.f.s were arbitrarily taken to be standard uniform. As one can see, the em- 
pirical levels of the tests based on S n and K n are reasonably close to the 5% significance 
level, except sometimes for strongly dependent samples (r = 0.75) in which case the tests 
are, overall, too conservative. The test based on S n displays worryingly high empirical 
levels for samples generated from the Clayton copula when r > 0.5 (in particular for 
d = 2), although, as expected from the theoretical results of the previous section, the 
situation improves as n increases. It is on the contrary too conservative, overall, for 
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Table 2: Percentage of rejection of Hq computed from 1000 samples of size n £ 
{50,100,200} generated under H 0jTn PI -ffi, c > where H\ c is defined in (20), k* = [nt\, 
C\ and C*2 are both either bivariate Clayton (CI), Gumbel-Hougaard (GH) or normal (N) 
copulas such that C\ has a Kendall's tau of 0.2 and Ci a Kendall's tau of r. 



CI GH N 



n 


T 


t 


Q 


C' 


qR 






C' 

>->n 


qR 




r< 


c 

>->n 


qR 




50 


0.4 


0.10 


7.5 


8.1 


5.7 


5.8 


6.1 


3.7 


4.3 


5.5 


6.1 


4.8 


4.3 


4.6 






0.25 


12.1 


10.6 


4.0 


8.4 


9.4 


5.1 


5.0 


8.8 


10.4 


7.6 


5.5 


7.7 






0.50 


18.0 


16.1 


6.3 


13.3 


12.1 


7.8 


4.8 


12.5 


12.3 


8.4 


5.8 


10.5 




0.6 


0.10 


14.5 


17.0 


5.4 


11.2 


11.4 


7.7 


6.2 


9.9 


11.4 


9.8 


7.2 


10.2 






0.25 


35.5 


34.4 


7.4 


28.0 


29.9 


21.4 


6.4 


27.5 


31.3 


21.4 


7.1 


29.2 






0.50 


47.3 


41.6 


7.0 


43.7 


45.3 


30.3 


8.9 


46.1 


46.0 


33.9 


9.1 


45.1 


100 


0.4 


0.10 


7.1 


8.7 


6.1 


6.7 


6.6 


5.1 


5.1 


7.0 


5.8 


5.2 


5.3 


6.0 






0.25 


18.8 


19.9 


5.2 


16.1 


16.9 


13.2 


5.8 


18.7 


14.9 


12.5 


6.1 


16.9 






0.50 


26.5 


26.4 


7.3 


26.7 


23.8 


18.8 


7.3 


27.4 


22.6 


19.1 


7.9 


24.8 




0.6 


0.10 


21.5 


25.1 


5.8 


19.6 


16.7 


12.0 


5.1 


18.3 


17.5 


16.9 


6.1 


19.1 






0.25 


65.1 


66.0 


6.3 


63.1 


61.2 


51.5 


7.5 


64.7 


62.9 


54.8 


9.9 


67.4 






0.50 


82.1 


81.6 


14.7 


83.8 


78.8 


69.7 


14.9 


83.3 


79.1 


73.7 


11.8 


86.0 


200 


0.4 


0.10 


11.1 


13.8 


5.8 


11.1 


8.3 


8.1 


5.5 


9.2 


9.5 


9.8 


5.0 


11.7 






0.25 


30.8 


33.9 


5.9 


31.5 


27.6 


24.8 


6.4 


31.6 


29.6 


28.3 


6.8 


34.8 






0.50 


47.1 


48.6 


9.0 


51.3 


45.8 


41.4 


8.7 


50.4 


47.1 


46.1 


9.4 


54.0 




0.6 


0.10 


36.4 


41.3 


6.8 


40.7 


34.4 


31.7 


7.1 


44.0 


36.0 


36.3 


6.7 


49.3 






0.25 


92.6 


93.2 


12.3 


93.8 


91.4 


88.9 


16.7 


95.0 


91.3 


90.2 


12.0 


95.5 






0.50 


98.9 


99.3 


22.2 


99.5 


98.5 


98.1 


22.0 


99.6 


99.3 


99.1 


21.1 


100.0 



samples generated from the Gumbel-Hougaard copula. The test based on S% appears to 
hold its level reasonably well in all the scenarios under consideration. 

To investigate the influence of negative dependence, a similar setting was considered 
in dimension two for normal and Frank copulas whose parameters were fixed so that the 
copulas have a Kendall's tau of r £ {—0.25, —0.5, —0.75}. The results (not reported) are 
similar to those obtained in the case of positive dependence. Finally, the empirical levels 
of the tests were computed from bivariate datasets of size n £ {50, 100, 200} generated 



from asymmetric copulas obtained as in (33) when C is the bivariate normal copula with 



parameter 0.97, a £ {0.25,0.5} and (3 = 1 (see also Figure [I]). The resulting rejection 
percentages (not reported) appear to be reasonably close to the 5% significance level in 
all scenarios. 



4.3 Rejection percentages under f/o,m H H\ c 

The aim of the first experiment under H m D i?i )C was to investigate the influence of 
a change in the strength of the dependence only on the power of the tests. Table [2] 
reports the rejection percentages of H computed from 1000 random samples generated 
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Table 3: Percentage of rejection of Hq computed from 1000 samples of size n e 
{50,100,200} generated under ifo.m H Hi, c , where H\ c is defined in (20), k* = [nt\, 
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under H Qrn fl Hi c , where H m is defined in ([6]) and ifi iC is defined in (20), k* = \nt\ for 
t G {0.1,0.25,0.5}, and C\ and C*2 are both either bivariate Clayton, Gumbel-Hougaard 
or normal copulas such that C\ has a Kendall's tau of 0.2 and C 2 a Kendall's tau of r G 
{0.4, 0.6}. As one can see, the tests based on S n , S n and n n consistently outperform the 
test based on whose power is relatively weak for the sample sizes under consideration. 
The test based on S n appears overall slightly more powerful than its two competitors 
for n = 50, while, for n = 200, the highest rejection rates are obtained, overall, for the 
test based on n n . Unsurprisingly, the estimated powers are the highest when the change 
occurs in the middle of the sequence. The conclusions remain very similar in dimension 
three (results not reported) except that all the rejection rates are higher. This could have 
been expected as an increase in dimension can be regarded as an increase in sample size 
for data generated from exchangeable copulas. 



As the theoretical hints of Section |3.2| were provided for copulas with positive depen- 
dence, it seemed important to empirically investigate the case of negative dependence. 
Table [3] is the analogue of Table [2] in that situation. This time C\ and C 2 are both 
either bivariate normal or Frank copulas such that C\ has a Kendall's tau of -0.2 and 
C 2 a Kendall's tau of r 6 {—0.4, —0.6}. While the rejection percentages of the tests 
based on S n , S n and n n remain comparable with those of Table we see that, somehow 
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Table 4: Percentage of rejection of H computed from 1000 samples of size n 6 {100, 200} 
generated under H 0tTn D Hi )C , where H\ c is defined in (20), k* = [nt\, G\ (resp. C2) is 
a <i-dimesional Clayton (resp. Gumbel-Hougaard) copula whose bivariate margins have a 
Kendall's tau of r. 
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surprisingly, those of the test based on S R are higher, though still substantially smaller 
than those of the three other tests. 

The next experiment under i?o,m H H\ c consisted of taking C\ and C2 from dif- 
ferent copula families but such that their bivariate margins have the same value r G 
{0.25, 0.5, 0.75} of Kendall's tau. In other words, the aim was to investigate the influence 
of a change in the shape of the dependence on the power of the tests. Because the test 
based on n n only measures a difference in dependence strength, it has no power against 
this type of alternatives. Note that this is a hard problem, as distinguishing copulas on 
the basis of low amounts of data is known to be difficult, for instance, from the corre- 
sponding goodness-of-fit testing literature (see e.g. Genest et al. , 2009 Kojadinovic and 



Yan, 2011). The fact that the change-point is unknown makes the problem under con- 



sideration even more difficult. Two of the most easily distinguishable copulas are the 
Clayton and the Gumbel-Hougaard because the former has lower tail dependence while 
the latter has upper tail dependence (see e.g. McNeil et al. , 2005, Chapter 5). The rejec- 



tion percentages for this scenario are reported in Table |4j As one can see, for the sample 
sizes under consideration, the test based on S R has hardly any power against this type 
of alternative. For the tests based on S n and S n , the highest rejection percentages are 
obtained for r = 0.75. Overall, the tests based on S n and S n appear to have roughly 
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Table 5: Percentage of rejection of Hq computed from 1000 samples of size n e 



{50,100,200} generated under H 0i1n n H\ >c , where H\ c is defined in (20), fc* = [nt\ 



and C\ and Ci are asymmetric copulas defined as in (33) from the normal copula with 



parameter 0.97. For d, (a,/3) = (1,6) and for C 2 , (a, 0) = (5, 1), with 5 G {0.25,0.5}. 
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equivalent power for r = 0.25. The fact that the rejection rates of the test based on 
S n are sometimes higher than those of the test based on S n when r > 0.5 should be 
interpreted with care as the former was found to be too liberal in the case of strongly 
dependent data in the previous subsection. Notice also that the results are not symmetric 
with respect to t — 0.5. Indeed, from the results for t = 0.25 and t = 0.75, we see that it 
seems easier to reject H when the largest of the two subsamples comes from the Clayton 
copula. Finally, we see that the estimated rejection rates tend overall to increase as d 
changes from 2 to 3, which could be explained, as earlier, from the exchangeability of the 
copulas. 

Similar experiments were carried for C\ and C 2 in other families. When one of the 
copulas was taken from the Clayton family, the tests based on S n and S n displayed some 
power for r = 0.75. In other situations, such as when C\ is from the normal family and 
C2 is from the Frank family, the tests had hardly any power for the sample sizes under 
consideration. Again, this is not surprising from empirical results obtained in the context 
of goodness-of-fit testing. 

As a last experiment under i?o,m H #i, C) we considered the situation in which C\ and 



C2 are asymmetric copulas defined as in (33) from the normal copula with parameter 



0.97. For Ci, we took (a,/3) = (1,5) and for C 2 , (a,P) = (5,1), with 5 e {0.25,0.5}, so 
that C\ and C2 have the same value of Kendall's tau but "opposite" asymmetry. The left 
(resp. right) plot of Figure [l] shows 10,000 realizations of C 2 for 5 = 0.25 (resp. 5 = 0.5). 
The results of the experiment are reported in Table [5j As expected, the power of the 
tests based on S n and S n are the largest when the asymmetry is the strongest (5 = 0.5). 
Again, for the sample sizes under consideration, the test based on has hardly any 
power against this type of alternative. 



4.4 Rejection percentages under Hq jC Pi H\^ m 

We concluded our Monte Carlo experiments by an investigation of the power of the tests 
for distributional changes solely due to a break in the first marginal c.d.f. Samples of size 
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n G {50, 100, 200} were generated under H Q c D Hi m , where H^ m is defined in (32), such 
that the k\ = \ntx\ first observations of each sample were from a <i-dimensional c.d.f. 
with copula C and N(0, 1) margins, and the n — \nt\\ last observations were from a d- 
dimensional c.d.f. with copula C whose first margin is N(/a, 1) and whose d-l remaining 
margins are N(0, 1). The values 0.25 and 0.5 (resp. 0.5, 1 and 2) were considered for 
t\ (resp. jj). The rejection percentages of H when C is a normal copula are partially 
reported in Table |6j Similar results were obtained for the other copula families considered 
in the simulations. As one can see, the test based on S„ is this time substa ntia lly more 



3.2 



because 



powerful than those based on S n , S n and K n . This was expected in Section 
of (28), which shows that the test based on S% is actually a general-purpose test for 
detecting distributional changes. As far as the tests based on S n , S n and K n are concerned, 
their overall weak power against this type of alternatives can be explained by the fact 
that they were designed for detecting changes in the copula. Their rejection rates are 
close to the 5% significance level for \i = 0.5 and are the highest for = (2,0.25) 

and t = 0.5. Notice that another reading of the results for \i = 0.5 is that the tests based 
on S n , S n and K n , regarded as procedures for testing directly H c , seem relatively robust 
against small breaks in mean in one margin. 



5 Concluding remarks 



The weak limit of the sequential empirical copula process C ra defined in ^ was obtained 
and the asymptotic validity of a multiplier bootstrap was established under the nonre- 



strictive smoothness conditions considered by Segers (2012). These theoretical results 



were then used to derive tests for change-point detection particularly sensitive to changes 
in the copula. Among several tests derived from the process D n defined in (21), the 



one based on the maximally selected Cramer- von Mises statistic S n defined in (24) was 



found to be the most powerful in finite-samples. It was compared with the analogue test 



based on the statistic S% defined in (27) derived from the process defined in (26) 



considered in Remillard (2010, Section 5.2) and Biicher and Ruppert (2013 Section 3.2). 
Theoretical hints suggesting that, at least in some scenarios, tests based on D n should be 
more powerful for detecting changes in the copula than tests based on were provided, 
and this was confirmed by extensive Monte Carlo experiments. In these experiments, 
the test based on S n and the one based on K n defined in (31) displayed, overall, the 
best finite-sample behavior. The latter was found, overall, to be the most powerful for 
detecting changes in the dependence strength. The former has however the advantage of 
being consistent against all alternatives under Hq j7U D Hi iC , where H 0jTn is defined in 
and H\ c is defined in (20). 



It is our understanding that none of the existing nonparametric change-point tests 
designed to be particularly sensitive to changes in the copula can be used to directly test 
Hq c defined in (JTl) unless H Q m defined in ^ is assumed to hold. If such tests are applied 
and H defined in ^ is rejected, a pragmatic approach would then be to additionally 
perform d univariate change-point tests on each of the component sequences. Should 
there be no evidence against i?o,m, rejection of H 0tC could be concluded. 

The latter approach is however no completely satisfactory and a future research di- 
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rection would be to design nonparametric change-point tests that can be used to test 



H 0;C whether H 0jTn holds or not (see also the discussion in Quessy et al. 2013, Section 



4). Another important future contribution would be to extend the approach considered 
in this work to the case of serially dependent observations. For that purpose, a good 



starting point is the work of Bucher and Ruppert (2013) 
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A Proofs of Propositions [T] and [2 

Proof of Proposition [T| To show the desired result, let us first show that 

d 



Ar, 



sup sup 



C n (s,u) 



Z n (s, u) 



0. 



se[o,i] u£[o,i] d 

It can be verified that A n can be rewritten as A n = maxi<fc<„ \Jk/nBk, where 



(34) 



B k = sup 

i*e[o,i] d 



c k (i,u) - \ z fc (i,«) - £)c^(«)z fc (i, «<>'>; 



(35) 



We can further write A n = max{ a/1/ nB\ , ^/2/nB 2 , max 3 <fc< n \Jk/nBk\. To show that 
A n — — > 0, it therefore suffices to show that max 3 <fc< n sJk/nBk 0. 

Let a n = n~ 1 / 4 (logn) 1 / 2 (loglogn) 1 / 4 , n > 3. Then, 



k n k B k 

max \ —Bk < max \ — a k x max — . 

3<fe<n V n 3<k<n \ n 3<k<n dfc 



First, notice that 



max \ —dk = n 

3<k<n V n 



-1/2 



max k 1/4 (log £;) 1/2 (log log A;) 1/4 = a n -)• 0. 



3<fc<n 



Next, (12), which is equivalent to lim sup B n /a n < oo with probability one, implies 



that maxs<k<n Bk/ak < sup A . >3 Bk/dk < oo almost surely, from which we obtain that 
max 3 <K„ \Jk/nBk 0. The desired result is finally a consequence of the weak conver- 
gence of Z n to Zc in £°°([0, l] d+1 ) and the continuous mapping theorem. ■ 



Proo f of P roposition [2| The proof starts as the proof of Lemma 2 of Kojadinovic| 



et al. 



(2011). Without loss of generality, we deal with For any u G [0, l] d , let w_i 



denote the vector (w 2 , . . . ,Ud) of [0, l] d 1 . Also, let 5 > be a real number such that 
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0<5<a<b<l — 5 < 1, and let n be sufficiently large such that, for any x G [a, b], 
x±rT x l 2 G [5,1 — 5]. It follows that, for any u\ G [a, b], uf n = Ui+n' 1 ^ 2 , u\ n = U\ — n -1 / 2 , 
and u~t n — itj~ n = 2n~ 1 / 2 . Also, for any u G [0, l] d such that U\ G [a, 6], we have that 



sup \C^\u)-C [1 \u)\ < sup 



uG[0,l] d 
ui G [a ,b] 



u£[0,l] d 
ui € [a, b] 



^ {C« B ,u-i) - Cfe,^!)} - <7N(«) 



2n-V2 

+ sup |C„(l,wJ n ,it_a) -^(l.tt^u-i) 

ue[o,i] d 



(36) 



where C„(l, ■) = y/n(C n — C); see also (pi). Since C' 1,1 ' exists and is continuous on 
V\ = {u G [0, l] d : < u\ < 1} according to Condition (C2), for any u G [a, 6] x [0, l] d_1 , 
there exists G (ux,uf n ) C [5, 1 — 5] and w7£ G (% n ,«x) C [5, 1 — 5], such that 

C«„,u_i) = CH+C [11 H«- 1/2 + C [1 ' 11 « n t ,i t _ 1 )(2n)- 1 

and 

C(%,„,«-i) = <7(u) - C7N(u)"~ 1/a + C^'VlS, "-i)^)" 1 , 
which implies that 



sup 

ixS[0,l] d 
ill £ [a, fa] 



2n- 1 /2 



^n- 1 / 2 sup |C^(u)| =0(n- 1 / 2 ), 



since C' 1 ' 1 ] is continuous on V - / = {it G Vi : U\ G [5, 1 — 5]} C Vi. 

For the second supremum on the right of (36), we start from (12) and write 

d 

C n (l,u) = Z n (l, «) - J] C b1 («)^n(l, + ^n(w), 

where sup u6 r 01 i<i \R n (u)\ = O (tit 1 / 4 (log n) 1 / 2 (log log n) 1 / 4 ) almost surely. Then, 
sup |C n (l,M^ n ,U_i) - C n (l,u^ n , U_i)| 

ug[0,l] d 



uiG[a,b] 



where 



< sup |C n (l,u) -C n (l,v)| < Ai >n + A 2jri + A 3)n + A 4 , 5 

I 11 !-"! I<2n -1 / 2 



A lin = sup |Z n (l, it) - Z n (l,v)| , 



|u 1 -i; 1 |<2n- 1 / 2 

A 2 , n = sup |CM(u)Z n (l,u«) - C [11 («)Z n (l,!;( 

| t i 1 -u 1 |<2n -1 / 2 

A 3 , n = ^ sup \{C®(u) - CM(v)}Z n (l,u^)\ 

j=2 u > v £Vf > u — i — w — 1 



24 



and A± >n = 2 sup ug [ ^ \R n (u)\ = O (n _1 / 4 (logri) 1 / 2 (loglogri) 1 / 4 ) almost surely. Pro- 
ceeding as in the proof of Proposition 4.2 of Segers (2012) for the term I n , it can be 
verified that, with probability one, 

A lin < sup \Z n {l,u)-Z n {l,v)\ =0(n- 1 / 4 (logn) 1 / 2 ). 

u.,vG[Q,l] d 
Vj ly|<2n -1 /2 

For A 2 - n , we have A 2 , n < B ljH + _B 2 ,n, where 

B hn = sup \C [1] (u) -C [1] {v)\ x sup |Z n (l,w (1) )| 

|u 1 - t ) 1 |<2n~ 1 / 2 



and 

^2,n = sup |C [1] (w)| x Ai )n = 0(n~ 1/4 (logn)~ 1/2 ) almost surely. 
For Bi n , from the mean value theorem, we obtain that, with probability one, 



B hn < sup \C™ (it)| X 



SUp^ ,l]^n(l^^)| x2n . 1/2(loglogn)1/2 = 0(n -l/ 2(loglogn) l/ 2)j 

(log logn) 1 ^ 



since sup n>3 sup ug [ |Z n (l, w^)) | /(log log n) 1 / 2 < oo almost surely by the law of the 
iterated logarithm for empirical c.d.f.s (see e.g. Chung, 1949, Theorem 2*). For the term 
^3,n, we proceed as in the proof of Proposition 4.2 of Segers (2012) for the term III n . 
Write A 3tTl = ^ =2 ^3,nj, where 



.4 



sup \{C [3] {u) -C [j] {v)}Z n {l,u U) )\ 



Tt,vGV^ > u — 1 — v — l 
|Ui —V-, ]<2n -1 / 2 



Let 7„ 



n 



- 1 / 2 (logn)(loglogn)- 1 / 2 , n > 3, let V> = {u E Vj : u 3 G [ 7n , 1 - 7 n]} C 



V^, and let V? n = [0,1] \ V? n . Then, using the fact that V{ is the disjoint union of 
Vf n V] n and 7/ n V^ 7 ", and that sup Uj „ e[0>1]d \ C [j] {u) - C [j] {v)\ < 1, we obtain that 



.4. 



3,n,j 



< A' 3tnd + Al nJ , where A', 



A" 



3,n,j 



su Puev ln |^n(l, u^) | and 



sup 

j 

u 1 -i> 1 |<2n -1 / 2 



«,»6V 1 s nl' 1 " ,u_x=w_l 



| - C w (w)} Z„(l, t* y 5) | . 



As ve rified in Segers ( |2012 , Eq. (4.3)) using Theorem 2 (iii) of Einmahl and Mason 
( |1988| ), A' 3 n j = O (n 1/4 (log n) 1/2 (log log ra) 1/4 ) almost surely. Also, from |Segers| ( |2012 [ 
Lemma 4.3), there exists a constant a > such that 

4Ui<« sup { '7"" 1 ' |z n (i,^))||. 

\<2n~ 1 / 2 

Using the fact that {uj( l-Uj)}~ l/2 < {7 n (l-7 n )} _1/2 = Q(7n 1/2 ) for all uj g [ y n , l-j n ], 
and that, as verified in [Segers (2012, page 775) for instance from [Mason (1981, Theorem 
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2), sup u6 r 01 ]d Z n (l, u^>) I {Uj(l — Uj)} x l 2 = O ((log n) 1 / 2 log log n) almost surely, we 
have, with probability one, 

A» >njj <a2n- 1 /*x sup J^'^lw X HXl ~ ln)}' l/2 = O (n-^(log\ogn)^) . 
ue[o,i] d i u A l - u i)} 1 

It follows that A 3>n = O (n~ 1//4 (logn) 1 / 2 (loglogn) 1 / 4 ) almost surely, which implies the 
desired result. ■ 



B Proof of Theorem U 



Let us first define additional notation. Let G n be the empirical c.d.f. computed from the 
unobservable random sample Ui, . . . , U n and let G n j, j G {1, . . . , d}, be the univariate 
margins of G n . Recall that, for any j G {l,...,d}, Rij >n is the rank of among 
Xy, . . . , X n j, or, equivalently, the rank of Uij among Uij, . . . , U n j. We also have that, for 
any i G {1, . . . , n} and j G {1, . . . , d}, G n j(Uij) = R^^jn = Uij >n (n + l)/n. Furthermore, 
for any j G {1, . . . , d}, let 

G-){u) = ini{v G [0, 1] : G n>j (v) >u}, u G [0, 1], 

be the quantile function associated with G n j. It is well-known that 



i i t/" fc:riij if (A; - l)/n < m < k/n, 

^ W 10 ifw = 0, 



(37) 



where Ui :n j < • • • < U n:n j are the order statistics associated with the jth coordinate 
sample Uij, . . . , U n j- Next, we define 



G- 1 (i») = (G- 1 1 («) l ...,G;i(«)) > ue[0,l]< 
Finally, for any m G {1, . . . , M} and (s, u) G [0, l] d+1 , let 



(38) 



s,tt) 



(m) 



i=l 



and 



2(m) ( 



|nsj 



[1 {Ui < - G n {G-\u)}] . 



(39) 



(40) 



i=i 



The proof of Theorem [T] is based on six lemmas. 
Lemma 1. Assume that (Ail ) holds. Then, for any m G {1, . . . , M}, 



sup sup 

se[o,i] ue[o,i] d 



4o, 



where zi m ' ) and zl m ' ) are defined in |i5|) and |55|), respectively 
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Proof. We have sup sg r 01 i sup u6 r 0)1 id Z„ (s, u) — Z { ™'(s, u 



j{m) 



< J n + K n , where 



J n = sup sup 

se[o,i] ue[o,i] d 



|nsj 



* i=i 



and 



K n = sup 
se[o,i] 



[nsj 



sup 

iie[o,i] d 



C M (u) -G LnsJ \ G^ s} (u) 



Now, 



and 



K n = max 

Kk<n 



In ^— ' 

i=l 



sup ICfeH-ClG, 1 ^)}! 
ue[o,i] d 



A; 

ICfcH-CjC, 1 ^)}! < ^|l(t/ i)fc <w)-l{C7i <G?(u)} , ue[0,l] 



i=l 



From (37), we have that, for any j 6 {1, . . . , d} and u G [0, 1], Uij < G k l -{u) is equivalent 
to (Rij,k — 1)A < M - Since Rij^/ik + 1) — (Rij k — > 0> ^ follows that, for any 
« G [0, \} d , 



d 



Rij,k 1 



Afrfc 1 ) . (4i) 

fc J fc+1 ' 



Consequently, 



sup |C fc (w) 
ue[o,i] d 



^{G^)}^ sup 



Rij,k 1 Rij,k \ , ^ 



(42) 



Hence, 



< — ^= max 

'Tl l<k<n 



1 ^ 



i=l 



(m) 



d 

< — = sup 

V n k>l 



(m) 



since, from the strong law of large numbers, sup fc>1 



1 k 

i=i 

fc _i Etid m) 



< oo almost surely. 
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It remains to show that J n — > 0. Using (41), we obtain that 

k 

i<k<n ue[0 y ]d y/n 

p 1 



J n <max sup ^^kH\HU iik <u)-l{Ui<G];\u)} 



< max sup — == 1 ( 

i<fe<n u6[0!l]d y/n ^ I 1^ V 

1 d 

i=i 



Rij k 
k + 1 



< — = > max sup > £ 



(m) 



n ^ i<fc<™ Uje[ o,i] . : 



fc + 1 



< —= y max max £ 

'71 ^— ' l<fc<n l<i<fc " 

- 1 /2 Tri fl.v^^_ lA (m) 



(m) 



< — p= max 

n l<fc<n 



^(m) 



The expectation of n x ' 2 max!< fc < n |^ mJ | is known to converge to zero prov ided E {(^ m) ) 2 } < 



Kosorok 



oo, which is the case since J OO {P(|£i" 1 ' > | > x)} 1 / 2 ds < oo (see e.g. 
10.5.1 and 10.5.2). It follows that J n — > 0, which completes the proof. 

Lemma 2. Assume that (A42) holds. Then, for any to 6 {1, ... , M}, 

Z^(s,u)-Z^\s,u) 



2008 



Exercises 



sup sup 

se[o,i] ue[o,i] d 



4o, 



where Z^ and Z^ are defined in (14) and (15), respectively. 



Proof. We have sup s6f0 xl sup 



*e[o,i] bu Pug[o,i]° 



zi m) ( s , w )-zi m) ( s , w ) 



< J n + K n , where 



J n = sup sup 

se[o,i] ue[o,i] d 



and 



K n = sup sup 

se[o,i] ■ue[o,i] d 



z( m )( s , w )-z( m )( s , w ) 



z( m )( s , w )-z( m )(,,^ 



We know that X„ — » from Lemma 



It remains to show that J n — > 0. 



It can be verified that J n < J' n + J", where 



j;=sup sup zH^uj-ZW^G-^n)} 
se[o,i] ue[o,i] d k - 1 



and 



4' = SU P 

*e[o,i] 



|ns] 

1=1 



sup 

ue[o,i] d 



{ <?£](«) \-C \G^ si (u) 



Now, 



J" < max 

Kk<n 



k 

Z^^ 



i=l 



sup \G k {u) -C{u)\ 

u<=[0,l] d 
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Using the law of the iterated logarithm for i.i.d. mean variance 1 sequences and the 



law of the iterated logarithm for multivariate c.d.f.s (see e.g. Kiefer, 1961; Dudley and 



Philipp , 1983), J" converges to zero almost surely since 



max 

3<fc<n 



sup \G k {u) - C(u)\ 
tte[o,i] d 



sup ue[0}1] d\G k (u)-C(u)\ l 
< max jf- — —pr x max , U9n — : ^ x —j= max log log k 

3<fc<n (fc log log k) 1 ' 2 3<k<n k^ 1 ' 2 (log log k) 1/2 y/Tl 3<k<n 

converges to zero almost surely. 
Concerning J' n , we have 

J' n = max sup \zl m \k/n, u) - Z^ {k/n, Gf(u)}\ . 



i*e[o,i] d 



Let a n = n 1,/2 (loglogn) 1//2 , n > 3. From the law of the iterated logarithm for univariate 

(43) 



c.d.f.s (see e.g. Chung, 1949, Theorem 2*), we have that, for any j e {1, . . . , d}, 

almost surely. 



lim sup 



sup ne[0 ,i] \G n ,j{u) -u\ l 



By considering the graphs of G n j and G n j, it immediately appears, as in 
Welner| (|1986[ page 86), that 



Shorack and 



sup \G n j(u) — u\ = sup \G n Au) — u\, 
ue[o,i] ' «e[o,i] 

from which we obtain that, for almost every u, there exists k w such k > k^ implies 

max sup \GIAu)(lo) — u\ < a k . 

l<j<d ue [ 0}1 ] 



(44) 



Hence, for almost every u, 



J' n (u) < max Y nk {u) + max Y nk {u) 

<k<n 



where Y n>k = sup 



ue[o,i]° 



Zr (k/n, u) - Zr { k/n, G k \u) } 

k u -l 



Now, 



max Y nk (u) < 
Kk<ku \/n 



u)\ -> 0. 



i=i 



For > we have 



Y n>k (u>) < sup |z£")(&/™,u)(u;) -Z^)(A;/n,«)(a;)| . 

■u,u6[0,l] d 

In other words, with probability one, J' n is smaller than 

L n =max sup \l^\k/n,u) -I^\k/n,v)\ 

l<k<n UtV£[0A] d 

Vj \v,j-Vj\<a k 
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plus a term that converges to zero almost surely. 

p 

To conclude the proof, it remains to show that L n — > 0. Let 



M n (a) 



sup 

•u., -!>£[(), l] d 

Vj I u j — Vj | < a 



1 n n 
i=l V i=l 



a > 0, 



where ^(«, «) = 1(17; < «) - C(v) - 1(17* < «) + C{v) and, for any / : [0, 1] M -> 
|o = sup UlWe [ ,i]d v)\. Then, L n can be rewritten as 



Vj I -u j — v a | < a 



L n = max 

Kk<n 



E^ m) Zj = max \ -. 
v i<fc<n V n 

1=1 a k 



max \l -M k (a k ). 

l<k<n V n 



Let v > 1. Then, 



E = E 



max 

l<k<n \ n 



v/2 



M k (a k ] 



< E 



£ 

,fc=i 



i//2 



< n max 

Kfe<n V n 



M k (a k 



E 



M fe (a fe ) 



From Lemma 2.3.6 of van der Vaart and Wellner (2000) with $(se) = x", x > 0, we obtain 
that, for any k > 1, 



E 



E 



"A- 



< 2^E 



1 



OA- 



It follows that 



where 



E(iZ) < 2"ra max 

l<fc<n \ n 



v/2 



MJa) 



n 

/n ' 
i=i 



E[{M fc (a fc )f] 



a > 0, 



is the oscillation modulus studied by Einmahl (1987). From his Theorem 5.2, we have 
that, for sufficiently large n, 

E[{M n (a)y} = 0{(a\oga- 1 y/ 2 ) 

when C is the independence copula. The proof of that result is based on Inequality 5.3 on 
page 73 of Einmahl (1987). As explained in Segers (2012, Appendix A), this inequality 
continues to hold for any copula C. Consequently, we have that 

E[{M n (a n )} u } = 0(b n ), 

where b n = n~^/ 4 (loglogn) !/ / 4 (logn) Iy / 2 . 
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Then, we write 



Tn max I -V^E \{M k (a k )Y] < 2 v n max (-\'\ k X max E^K)}' 

3<fc<n \n J 3<k<n \n J 3<k<n b k 

The last maximum is bounded, while 



v/2 



/k\ U/2 

2 v n max - b k = 2 u n l - //2 max k v/2 b k = 2 u n 1 - u/4 (loglognY /A (lognY /2 , 

3<fc<n / 3<A:<n 

which converges to zero for v > 4. It follows that E (LJQ converges to zero for any v > 4, 

p 

which implies that L n — > 0. ■ 
Lemma 3. Assume that (A42) holds. Then, for any m G {1, . . . , M} and j G {!,..., c?} ; 



sup sup 

se[o,i] tte[o,i] d 



Ao, 



where CuL defined analogously to (11), and Z^ mj zs defined in (14)- 
Proof. For any 5 G (0, 1/2), we have 



7 (m) 



sup sup 

se[o,i] tie[o,i] d 



where 



J n = max 

Kk<n 



k 1/4 

n l/8 



sup 

u6[0,l] d 
jjG[«,l-«] 



C7W(ti)-CW(tt) 



" ue[o,i] d 



zi m) (i,«W) 



and 



K n = 5 sup sup \Z^(s,u^)\ 

se[0,l] «e[o,i] d 



since Cn (tt) < 5 for all w G [0, l] d and n > 1 as shown in Kojadinovic et al. (2011 
Proof of Proposition 2). Now, since, for any s G [0, 1], zi m ^(s,tt) = when u is on the 
boundary of [0, l] d , 

sup sup |Z^(s,u^)| < sup sup \Z^\s,u^) -Z^\s,v^)\ 

se[0,l] ™e[o,i] d se[0,l] «,»6[o,i]' i 

Uj£[S,l — S] \uj—Vj\<5 

< SUp \% { n n \s, U^) - Z { ™\t, V®)\ , 

s,te[o,i],u,D6[o,i] d 



Is-H+piuj ,Vj)<V S+S 2 



where 



p 2 (u, v) = E [{1(U < u) - u - 1(U < v) + v} 2 ] = \u - v\ - (u - vf 
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and U is a standard uniform random variable. The weak convergence of implies the 
asymptotic equicontinuity of (s, Uj) >->■ Zn"^(s, 1, . . . , 1, Uj, 1, . . . , 1) with respect to the 
natural semimetric on [0, 1] x [0, 1] (see van der Vaart and Wellner, 2000, Chapter 2.12), 
i.e., for any e > 0, 



limlimsupP I sup {s , u ij) ) - Z { ™\t, v ij) ) \ > e \ = 0. 

74-0 n-toc | «,te[o,i],u,«e[o,i] d 

\»-t\+p(uj ,Vj)<y 

Hence, for fixed e > and r\ > 0, we can choose 5 > sufficiently small so that, for 
sufficiently large n, P(K n > e) < rj/2. 

Let us now show that, for the chosen S, J n converges in probability to zero, which will 
imply that, for sufficiently large n, P(J n > e) < i]/2 and will complete the proof. 



From the law of the iterated logarithm stated in (43), we have that, for any j G 
{!,...,<*}, 

sup u6[0 ,i] n 1 ' 2 \G nd {u)-u\ 



SU P n i m/2 

„> 3 (log log n) 1 ' 1 



< oo almost surely. 



Since \l{U n j < u) — u\ < 1 for all u G [0, 1] and n > 1, for any v > 0, 

{sup ue[0 ,i] |l(£/ ni < u)-u\} v 



E 



sup 

n>3 



( log \ogn) u l 2 



< oo, 



which, from Corollary A. 1.8 of van der Vaart and Wellner (|2000|), is equivalent to 

E 



sup 

n>3 



{sup ue[0jl] n l ' 2 \G n j{u) - u\Y 



(log logn) 1 '/ 2 



< oo, 



and implies that 

E [{ s up M6[0 ,i] n l l 2 \G n j{u) -u\} u ] 



sup 

n>3 



(log \ogn) v / 2 



< E 



sup 

n>3 



{suPng[o,i] n 1/2 \G nJ (u) - u\} 1 
(log \ogn) u / 2 



< oo. 



Let v > 1. Applying Lemma 2.3.6 of van der Vaart and Wellner (2000) with $(x) = x u , 
x > 0, we obtain 



sup E 

n>3 



su Pue[o,i] a 



(log log7?,) y / 2 



< sup E 

n>3 



(log logTi)^/ 2 



< OO. 

(45) 



Let us now show that J n — > 0. We have 



max 

3<fc<n 



k 1/4 



sup 

ie[o,i] d 



Cf{u)-C^\u) 



sup 



n3/8 ue[o,i] d 



< L n x M n , 
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where 



Cf{u)-C^{u) 



max sup - — — — . - j- - — — — r . x n 1 ^ 8 max (lo£r/c) 1 ^ 2 (lo£rlogfc) 3//4 



and 



M n = max 

3<k<n 



n tt6[0,l] d 



The first maximum in L n is strictly smaller than oo almost surely by Proposition [2j which 

as P 

implies that L n —4 0. It remains to show that M n — > 0. For v > 1, we have 

k^ 4 



E (M£) = E max 

3<fc<n 



' i ltte[o,i] d 



Z^(1,«W; 



< E 



£ 

,fc=3 



fc^ 4 



3i//8 



sup 

ue[o,i] ci 



zJ B) (i,«w; 



< n max - ,„ E 



3<fc<n n 



3v/8 



sup 

tte[o,i] d 



Z< m) (l,u<fl) 



< n 1 " 3 ^ 8 max r /4 (loglog fcW 2 x max 

3<fc<ra 3</c<n 



[{ 



E lisup u6[0jl]d 



(log logk) u / 2 



The second maximum on the right of the last inequality is bounded by (45), while the 

first maximum is equal to n 1_I/ / 8 (loglogn) I///2 , which converges to zero for v > 8. Hence, 
p 

M n — > 0, which completes the proof. ■ 
Lemma 4. Assume that (Ml ) holds. Then, for any m e {1, . . . , M}, 



sup sup 

se[o,i] ue[o,i] d 



z( m )( s ,n)-i( m )( s , w ) 



4o, 



where zi m ' ) and zi m ' ) are defined in (16) and (40), respectively. 

Proof. The proof is a simpler version of the proof of Lemma [T] and is omitted. 
Lemma 5. Assume that (Ml) holds. Then, for any m £ {1, . . . , M}, 

sup sup I Z^ m) (s, u) - Zl m) (s, u) | 4 0, 
se[o,i] ue[o,i] d 



where zi"^ and zi"^ are defined in (14) and (16), respectively. 
Proof. We have sup se[01] sup ue[0]1]d 



zi m) ( s , w )-zr j (s,u 



(m) 



< J n + fC, where 



J n = sup sup 

se[o,i] ue[o,i] d 
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and 



K n = sup sup 

se[o,i] ue[o,i] d 



We know that K n — > from Lemma 



It remains to show that J n — > 0. 



It can be verified that J n < J' + J" where 



j;=su P sup (z^^^-z^i^g; 1 ^)}! 

se[o,i] we[o,i] d 



and 



J n = sup 

se[a,i] 



[nsj 



In / — ' 

i=l 



x sup ^{G^Hl-ClG^H}!. 

iiS[0,l] d 



Now, 



J'n < SUP 

se[o,i] 



[ns\ 
i=l 



x sup \G n (u) — C(u)\ — y 
■ue[o,i] d 



since the first supremum converges weakly to the supremum of a Brownian motion and 
the second supremum converges almost surely to zero by the Glivenko-Cantelli lemma. 



Concerning J' n , we have, using (43), that, for n sufficiently large 



J' n < sup sup \Z ( ^\s,u)-Z^ n \s,v)\, 

se[0,l] u,v£[0,l] d 

Vj \ uj —Vj\<a n 

where a n = n™ 1 / 2 (log log n) 1 ^ 2 , n > 3. The quantity on the right converges in pr obability 
to zero due to the asymp totic equicontinuity of Indeed, as discussed in 



Vaart and Wellner 



stated in terms of 
5 n j. 0, as 



van der 



(2000J Chapter 2.12), the asymptotic equicontinuity of zi can be 
;he natural semimetric on [0, 1] x [0, for any decreasing sequence 

sup \Z^\s,u)-Z^\t,v)\ 4 0, 

S,tS[0,l],Ti,DS[0,l] d 

where p 2 c {u, v) = E [{1{U < u) — C{u) — 1(U < v) + C(v)} 2 ], and U is a random vec- 
tor with c.d.f. C. 

Now, it is easy to verify that 

p 2 c {u,v) = \C(u) — C{u Ad) + C{v) — C{u A v) - {C(u)-C(v)} 2 \ , 
which implies that 



p 2 c (u, v) < C{u Vd) - C(u Av) + {C(u) - C(v)f 



< d max \uj — vA + d 2 max \uj — vA 

l<j<d ±<j<d 



since C satisfies the Lipschitz condition 



\C(u)-C(v)\<J2\^-^\, u,ve[0,l} d . 



(46) 
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Hence, with 5 n = \J da n + d 2 a^, 



J' n < SUp SUp 

*6[0,1] v.,ve[o,i] d 

Vj \Uj —vA<a„ 



\Z^\s,u)-Z^\s,v)\ < sup |zM(s,u) -Z^\s,v)\ 

s€[0,l],ii,u6[0,l] d 
p c (v.,v)<5 n 

< sup \Z^ n \s,u)-Z^ n \t,v)\^0, 

s,t£io,i],u.,veio,i] d 

\s — t\+Pfj{-u.,-u)<S n 

which completes the proof. ■ 
Lemma 6. Assume that holds. Then, for any m G {1, . . . , M} and j G {1, . . . , d}, 

sup sup [\C®(u)-C®(u)\ |Z^(s,«^)|] 4 0, 
se[o,i] ■ue[o,i] d 



where Cn' and Z)^ 1 ' are defined in (TiTl) and (14\), respectively 



7 {m) 



Proof. For any 6 G (0, 1/2), we have 



sup sup [\C^(u)-C^(u)\\Z^(s,u^)\]<J n + K n 
se[o,i] ue[o,i] d 



where 



and 



J n 



sup \C®(u)-C®( 

ue[o,i] d 

« ,-6[«,l-«] 



uj| sup sup \zW(s,u®) 
se[o,i] tte[o,i] d 



K n = 5 sup sup |Z^(s, it 



se[o,i] iis[o,i] d 



since Cn\u) < 5 for all u G [0, l] d and n > 1. Proceeding as in the proof of Lemma|3j for 
fixed e > and rj > 0, we can choose 5 > sufficiently small so that, for sufficiently large 
n, P(K n > e) < r]/2. For that 5, J n converges in probability to zero by Proposition p] 



and since sup se[0>1])Ue[0il]<i 



Op(l). It follows that, for sufficiently large n, 



P(J n > e) < i]/2, which completes the proof. ■ 
Proof of Theorem Q3 We only prove (ii) as the proof of (i) is similar. A consequence 



of the sequential extension of the multiplier central limit theorem proved in |Holmes et al. 
fl2013l Theorem 1) is that 



(z m z«..,zW)-(z C) zi;)... ) 4 M) ) 



(47) 



in {£°°{[0, l] d+1 )} A/+1 , where Z n and zi m) are defined in (tal) and dl4fc, respectively, Z c is 



a C-Kiefer-Miiller process and z9,\ . . . , Z^ 1 ' are independent copies of Z, 



(M) 



Now, for any m G {!,..., M}, let 



C£ m >( a , u) = Z^\s, u) -Y^d®(u)I^\a,vP), (s,u) G [0, 1] 
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Then, from the continuous mapping theorem and the fact that A n defined in (34) con- 
verges in probability to zero, we have that 

(C n X£\.-.Xl M \c£\...,CW) - (c c ,C«... ! C( M » 1 C2 ) ,... 1 Clf ) ) (48) 

in {£°°([0, i]rf+i)}2M+i ) where £c is defined in and Cg } , . . . ,C[f ) are independent 
copies of Cc- Next, from Lemmas [3] and |6j we obtain that 

C n , . . . , . . . , e«) - (c c> C«, ... , Cjf \ cg\ . . . , C<f > 

in {£°°([0, l]^+i)}2M+i ) w here, for any m G {1, ... , M} and (s, u) G [0, 



3=1 



and 



,0> 



3=1 

The desired result finally follows from Lemmas [2] and [5J and the fact that C$ (it) < 5 for 
all u G [0, l] d and n > 1. ■ 



C Proofs of Propositions [3], [4] and [5] 

Proof of Proposition [3| Under Hq, D n can be rewritten as 

O n (s, «) = {1 - A n (s)} C n (s, u) - A n (s)C*(s, u), 
where C n is defined in rt3J) and where, for any (s,u) G [0, l] d+1 , 



C* n (s,u 



) = ^{i-A n ( S )}{c H w-cw} = ^ E {i(t/* n _ M < W )-c(u)} 

(49) 



i=[ns\ + l 



Let us first show that 



A* n = sup sup 

se[o,i] ite[o,i] d 



c*( a ,„)_-L £ < «) - c(u)} 



u 



i=l v ' ' i=\ns\+l 

Proceeding as in the proof of Proposition [TJ we have that 
A* n = max \jn — k/nB* 



4 0. (50) 



0<k<n~\ v n k 1<k<n v 
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where 
B 



1 



n-k ~ SU P / r 

ue[o,i] d yn- k 



n n 

E {!(^ < «) - c(«)} - E t 1 ^ ^ «) - 



i=fc+l 



+Ec b1 H x E W 17 * - uU) )- u i} 

3=1 



i=k+l 



k G {0,...,n- 1}. 



Now, clearly, for any fc G {1, . . . , n}, (U hk , U kjk ) and (U*_ k+1>k , U* k ) have the 



same distribution. As a consequence, for any k G {!,••■ , w}, -B£ and 5^ defined in (35) 



have the same distribution. More generally, (Bi,...,B n ) and (B*, . . . , B*) have the 
same distribution. Indeed, (B*, . . . , B*) is nothing else than the version of (£>i, . . . , B n ) 



computed from the sequence X n , X n _i, . . . , X\. It follows that A* and A n defined in (34) 



have the same distribution. Consequently, A n —4 implies that A* n — > 0. 

Next, from the weak convergence of Z n to Z c in e°°([0, l] d+1 ) and the continuous 
mapping theorem, we have that (s, u) i-> Z n (s, it) — Ylj=i C^(u)Z n (s, it^) and 



(s, u) i-> 



1 n d 1 n 

-= e {i(^<^)-^)}-E^' ] ( w )x^ E w^<« (, ' 3 )-«i} 

^ i=|nsj + l j=l * i=|rasj+l 



= Z n (l,«) - Z n (s,-u) - Ec b1 N {Z n (l,w a) ) -Z n (s,«^)} . 
jointly converge weakly to Cc and C^, respectively, in {£°°([0, l] d+1 )} 2 . The desired 

as P 

result is finally a consequence of the fact that A n —4 and A* — > 0, and the contin- 
uous mapping theorem. The expression of the covariance function follows from simple 
calculations. ■ 

Proof of Proposition [4| We only prove (ii) as the proof of (i) is similar. For any 
m G {1, . . . , M} and (s, u) G [0, let 



3=1 



where zi"^ is defined in (|14|). From (|47|), the continuous mapping theorem and the fact 



that A* n defined in (50) converges to zero in probability, we obtain that 



(51) 

in {£°°([0,l] d+1 )} 2M+1 , where C* is defined in @, Q(s,u) = C c (l,v) - C c (s,«) for 
all (s, it) G [0, and C£»^, . . . , C^ M ' are independent copies of C£. 

Next, let us show that, for any m G {1, . . . , M}, 



J n = sup sup 

se[o,i] ue[o,i] d 



lzi m) (i,w (i) ) -zi m) (s,w w : 



4o. 
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It is easy to verify that J n , written as a maximum over 1 < [ns\ < n, is nothing else than 



the version of sup sg [ ^ sup u6 [ 01 ] a 



C\ j] ,(u) 

Ins] V / 



Cb1( 



u 



computed from the 



sequence {^ n m \ X n ), X n _i), . . . , (C,[ m \ Xi), which implies that the latter quantity 

and J n have the same distribution. Lemma [3] therefore implies that J n A- 0. From the 
latter fact and Lemma |6j we therefore obtain that 



c*, b/ M \ t^\. . . , - ( c;, c;:" c;v -e;- c;, 



(n*,{M) ^,(1) fr*,(M) 
-C ' • • • ' ' ' • 

in {£°°([0, l]rf+i)}2M+i ; where, for any m G {1, ... , M} and (s, u) G [0, 



^*,(m) 



>, «) = ZM(1, «) - Z^s, u) - £ C^J nsJ («) {ZS»)(1, u«) - Zi m )( S , 



i=i 



and 



Now, reasoning as previously, for any m G {1, . . . , M}, 



= sup sup 

«e[o,i] ue[o,i] d 



ZM(l,u)-ZW(s,u)-Z£ (m) M 



is the version of sup sg [ 01 ] sup ue [ 0)1 ]d Zj^(s, it) — Z^ (.s. u) 



?( m ), 



computed from the sequence 



(£n , X n ), (Ci-i) ■X'n-i );•■■) (£i m \ which implies that the latter quantity and K" n 

have the same distribution. Lemma [2] therefore implies that K n — )■ 0. From the latter 
fact, Lemma y and since Cn\u) < 5 for all u G [0, l] d and n > 1, we thus obtain that 



*,(M) 



Q, Q« . . . , QW, QW, . . . , QM) - (C* c , C*/\ . . . , C£ w , C^, ...,C c 

in {£°°([0,l] d+1 )} 2M+1 . 

From (|47)), the continuous mapping theorem and the fact that A n and A* n (defined 



in (34) and (50), respectively) converge to zero in probability, we have that the weak 



convergences stated in (48) and (51) occur jointly. This implies that the weak convergence 



stated in the previous equation occurs jointly with that stated in Theorem [I] (ii). The 
desired result finally follows from the continuous mapping theorem. ■ 

Proof of Proposition [5| We first prove (i). Recall that U\, . . . , U n is the unobservable 
sample obtained by applying probability integral transforms to the coordinate sample of 
Xi, . . . , X n , and, for any (s, u) G [0, let 

B n (s,u) = V^Us) {1 - K(s)} [G LH {G^ sJ (u)} - G^ LfMj {G5:f H (tt)} 
where, for any k G {1, . . . , n — 1}, 



G,(u) 



1 k 



< u) 



and 



G 



n—k 



U 



1 n 

— r E !(^< u ^ [0,1]' 



i=l 



n — k 



i=k+l 
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with the convention that Gq(u) = and Gq(u) = 0, and where, for any k G {1, . . . , n— 1} 
and u G [0, l] d , G k \u) and G£ .fc(tt) are defined analogously to (p8h. 

Now, n~ 1/2 sup se[0il] sup ug[01]d \B n (s, u) - D n (s, u)\ < A n + B n , where 

k(n — k) 



A„ = max 



l<fc<n n 



sup \C k (u)-G k {G k \u)}\ 



u6[0,l]' 



and 



B n = max 



k(n — k) 



l<k<n n z 



sup \C*_ k (u) - G* n _ k {G^_ k (u)}\ . 
ue[o,i] d 



Then, from (42), 



k(n — k) d d kin — k) d d 
A n < max — < — and B n < max < 



i<k<n n 2 k n 



which implies that n 1//2 sup se r Q1 i sup ug r 01 ]d 
to show (i), we can show the result with D n replaced by 
Next, notice that, for any s G [0, 1] and u G [0, l] d , 



i<k<n n 2 n — k n 
(s, it) — 3D) n (s, it) | < 2d/n. It follows that, 



K t (s,u) 



s(l - s)Ci(ti) - s{(t - s)d(u) + (1 - t)C 2 (n)} if s < t, 
(1 - s){tCi(«) + (s- t)C 2 («)} - s(l - s)C 2 (u)} if s > t. 



Let us therefore split the supremum over s according to the cases s G [0, t] and s G [t, 1]. 
Then, we write sup sg[ot] sup ug[01 ]<i |n~ 1/2 D„(s, u) - K t (s,u)\ < I n + J n , where 



sup sup 

se[o,t] ue[o,i] d 



A n (s) {1 - X n (s)} G [nsi {G^ s} (u)} - s(l - S )d(u) 



and 

J n = sup sup 

s6[0,t] ue[0,l] d 

Now, 



A n (s) {1 - \ n (s)} G*_ M KC L Lj(«)} - - + (1 " 0C 2 (u)} 



/„ = sup sup 

se[o,t] ue[o,i] d 

< sup sup 
se[o,t] ite[o,i] d 



A n ( S ) {1 - X n (s)} [G lnsl {G^ si (u)}-C 1 (u)\ +o(l) 



+ max sup 

l<fc<|nt| ue [ ,i]d 



k(n — k) 



rr 



[Ct{G-j;\u)} - c x (ii 



where Z n is defined in (|9]). The first supremum on the right of the last inequality is smaller 
than sup sg [ i ] sup ug [ 0)1 ]d |n~ 1//2 Z n (s, it)| because G7 ,(it) G [0, l] d for all (s,it) G [0, l] d+1 . 
It therefore converges to zero in probability as sup sg r ot i sup ug [ 01 ]d |Z n (s, u)\ converges in 
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distribution. From (46) and (44), the second term on the right of the previous inequality 

su P« je [o,i] \ G k,j( u j) ~ u i\ 



restricted to k > 3 is smaller than 

d 



max — sup \G, — uA < max — a k x 



3<fc<KJ n ^—[ Uj eio,i} 



3<k<[nt\ Tl 



max 

3<fc<|ra*J 



a k 



< a n x d max sup 

l<j<d fc>3 



su P %e [o,ij l g fc,i( M j) ~ u j\ 



where a k = k 1//2 (loglog k) 1 ^ 2 , k > 3. By (43), the second term in the product on the 

right is strictly smaller than oo almost surely. Since the first term converges to zero, we 
p 

obtain that I n — >• 0. 

p 

Let us now show that J n — > 0. We have J n < L n + M„, where 



sup sup 

s£[0,t] ue[0,l] d 



An(a) {1 - A„00}C^ LH {Cfcf„j(u)} 

(t - S )^{G-Cf nsJ («)} + (1 - ^{^^(n)} 



and 

M n = sup sup 
se[o,t] ue[o,i] 



8 ((t - s) [ClKC^j («)} - Cl(tt)] + (1 - f) [^{G^j («)} " Cb(«)] ) 



Now, 



L n < sup sup |A n (s){l-A n (s)}G;_ LnsJ (w)-s[(t-s)C 1 (w) + (l-t)C 2 (w)]| 

se[0,i] -itG[0,l] d 



fc(n — fc) 

< max sup 

i<k<[nt\ n 2 ue[o,i] d 

< L n>1 + L n> 2 + o(l), 



G* n __ k (u) - - d(u) - L J C 2 (u) 



n — k 



n — k 



o(l) 



where 



k(\nt \ - k) 
L n l = max sup 



l<fc< \nt\ 



n 



ne[o,i] d 



[nt\ 



[nt\ - k 



i=k+l 



<n x l 2 sup sup \L n (t, u) — Z n (s, u) 
se[o,t] tte[o,i] d 



and 



L n ,2 = sup 

ii6[0,l] fl 



rrr J2 l(U,<u)-C 2 [xi, 

n — \nt\ 

L J i=[nt\+l 



x max 



k(n — \ nt\ ) 



\<k<[nt\ n z 



From the continuous mapping theorem, we obtain that L nji — > 0. Also, t/^j+i, . . . , U n 
being a random sample from C 2 , we have, from the Glivenko-Cantelli lemma, that 



L„2 — — > 0. Hence, L n — )■ 0. 
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Using (46) and (43) as previously, we obtain that M n — ^ 0, which implies that 
J n A 0, and therefore that sup sg [ 0jt ] sup ug r 01 id \n~ 1 / 2 'B n (s, u) — K t (s, u)\ A- 0. 

Proceeding similarly, one can show that sup sg [ t l ] sup^g^ ^d |n~ 1//2 D n (s, u)—K t (s, u)\ — > 

0, which implies that sup sg j 01 ] sup^g ^d |n~ 1 / 2 D n (s, u) — K t (s, u) \ — > and therefore the 
desired result. 

Let us now prove (ii). For any m G {1, . . . , M}, let 

\ns\ 



w (m) 



4= $^d m) {l(^ < «) - (a,u) G [0, 1] 



and 



Wl m) ( 



|n*j 

i=i 



From the proof of Theorem 3 of 



Holmes et al. 



(2013), sup sg[011 sup 



e[o,i] &u Pwe[o,i] 



d+l 



d W„ '{S.U 



Op(l) and sup sg [ 01 ] sup^g^^jd \W n (s,u)\ = Op(l) for all m G {1,...,M}. Then, from 
Lemma [TJ the fact that 



ZiT\s,u) =W^{s,G^ sj (u)}, (s,u) G [0,l] d+1 , 



where zi" 1 ^ is defined in (39), and the fact that 



sup sup 

se[o,i] ue[o,i] d 



< sup 


sup 




se[o,i] 


tte[o,i] d 





we obtain that sup sg [ 01 ] sup ug [ 01 ]d |Z« = Op(l). Hence, sup sg [ 01 ] sup Mg r 01 id {ClT'is 



(m). 



Op(l) since Cn\u)< 5 for all u G [0, l] d , j G {1, . . . , d} and n > 1. Similarly, using (40), 
wi™' and Lemma 



we obtain that sup sg [ 01 ] sup ug [ 01 ]d \Cn(s,u)\ = P (1 



Proceeding similarly, it can be verified that the corresponding results for C^ m ' ) and 
^qJ^^ which concludes the proof. ■ 
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